Extrinsic spin Hall effect in graphene 
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We show that extrinsic spin Hall effect can be engineered in monolayer graphene by decoration 
with small doses of adatoms, molecules or nano-particles originating local spin-orbit perturbations. 
The analysis of the single impurity scattering problem shows that intrinsic and Rashba spin-orbit 
local couplings enhance the spin Hall effect via skew scattering of charge carriers in the resonant 
regime. The solution of the transport equations for a random ensemble of spin-orbit impurities 
discloses that giant spin Hall currents are within the reach of current state-of-the-art in device fab- 
rication. The spin Hall effect is robust with respect to thermal fluctuations and disorder averaging. 
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The spin Hall effect (SHE) [T]^, that is, the appear- 
ance of a transverse spin current in a non-magnetic con- 
ductor by pure electrical control has been predicted to 
occur in materials with large spin-orbit interaction (SOI). 
Over the last decade, its study has lead to an intense ex- 
perimental activity [SHH] , due to its potential application 
in spintronics. Recently, the SHE has been explored for 
replacing ferromagnetic metals with spin injectors in ap- 
plications [ini E] , opening the doors to the development 
of spintronic devices without magnetic components. 

The activation and control of spin-polarized currents 
in solid-state is both of fundamental and technological 
interest. SHEs could be used for an efficient conversion 
of charge current into spin-polarized currents. The ratio 
of the spin Hall current to the steady-state charge cur- 
rent, commonly known as the spin Hall angle ^sh, mea- 
sures this efSciency and it is the most important figure of 
merit for practical applications. Broadly speaking, SHEs 
in metals and semiconductors originate from (i) extrinsic 
mechanisms, which are due to spin-dependent scattering 
of electrons by impurities in the presence of SOI [TH3], 
and (ii) intrinsic mechanisms, entirely due to SOI in the 
electronic band and that occur in the absence of any scat- 
tering process. In semiconductors, the spin Hall angles 
are in the range of 0.0001 - 0.001 Q. On the other 
hand, the experimental values of Ogn for metals can be 
considerable larger, being of the order of 0.01 for Pt [12] 
and 0.1 in a recent measurement performed in Ta [TT] . 

After its successful isolation, graphene [13] has also be- 
come the subject of intensive study in spintronics [I4l418j . 
In this material, electrons can propagate ballistically and 
carrier density and carrier polarity can be controlled by 
an external gate. Spin-orbit and hyperfine interactions 
are extremely weak in graphene, and therefore the spin 
coherence length is expected to be long [THl [2D] . These 
characteristics make graphene appealing for passive spin- 
tronic applications e.g., as a high-fidelity channel for spin- 
encoded information [3T]. One striking possibility is to 



manipulate electronic properties of graphene for active 
spintronics. This may be achieved via spin-orbit split- 
ting of the band dispersion e.g., by bringing heavy metal- 
lic atoms in close contact to graphene [22l, or by locally 
inducing sizable SOI (~ 10 meV) ^i24j; in Ref. [23], 
distortions induced by covalently bonded impurities were 
predicted to produce the desired effect, and Ref. [23] sug- 
gests local SOI enhancement via tunneling of electrons in 
and out a heavy atom. Moreover, it has been proposed 
that, in the presence of SOI, graphene could exhibit the 
quantum spin Hall effect [35]. 

In this Letter, we consider a monolayer of graphene 
decorated by a small density of impurities generating a 
spin-orbit interaction in their surroundings. We show 
that robust SHEs develop through asymmetric (skew) 
scattering events. Crucially, and unlike two-dimensional 
electron gases (2DEG), for which resonant enhancement 
of skew scattering [26) requires resorting to fine tuning 
and sometimes to exotic phenomena such as the Kondo 
effect [371 US] : our proposal takes advantage of graphene 
being a atomically thin membrane, whose local density 
of states easily resonates with several types of adatoms, 
molecules, or nano-particles. Resonant scatterers have 
been predicted to play an important role in charge trans- 
port at high electronic densities [5S1 |3D]. Here, we ar- 
gue that a similar physics is behind a huge potential of 
graphene for extrinsic SHEs. The decoration with small 
doses of adatoms/nano-particlcs only partially suppresses 
the charge carrier mobilities of graphene devices, which 
combined with large spin diffusion lengths and Fermi en- 
ergy tunability, makes this material a promising candi- 
date for spintronic integrated circuits with SHE-based 
spin-polarized current activation and control. 

According to our calculations, extrinsic spin Hall effect 
in graphene, as that recently reported in hydrogenated 
graphene [31j . can originate from skew scattering alone. 
The latter is absent in the first Born approximation [32] , 
and therefore we compute transport relaxation rates non- 




FIG. 1: Schematic picture of extrinsic spin Hall effect gener- 
ated by transport skewness. An impurity (sphere) near the 
graphene sheet causes a local spin-orbit field with range R. 
The scattering of components with positive (negative) angu- 
lar momentum is enhanced (suppressed) for charge carriers 
with Sz = 1 {sz ~ —i), resulting in a net spin Hall current. 



perturbatively via exact partial-waves expansions. We 
predict that clustering of adatoms or nano-particles leads 
to an enhancement of the spin Hall angle comparable to 
those found in pure metals {0sk ~ 0.01 - 0.1 [51 171 [12]). 
In order to study the amount of extrinsically gener- 
ated SHE and its dependence on impurity type, Fermi 
energy and temperature, we consider a continuum model 
of graphene decorated with a small concentration of lo- 
cal SOI impurities with different symmetries (see later). 



The starting point is the continuum-limit Hamiltonian of 
graphene Ho = Hvpirza^Px+c^yPy), where p = {px,Py) is 
the 2D kinematic momentum operator around one of the 
two inequivalent Dirac points {K and K'), vp « 10^ m/s 
is the Fermi velocity, a and r denote Pauli matrices, 
with (Jz = ±1 [tz = ±1] describing states on A(B) sub- 
lattice [at K{K')]. The spin-orbit splitting in the band 
structure of pristine graphene is of the order of 10 /xeV 
and therefore can be safely neglected [50]. The impu- 
rities play the role of scattering centers that locally in- 
duce SOI of the intrinsic type Vso = ^i{t^)tzO'zSz and/or 
Rashba-type Vso = ^RiT^)iTzCrxSy — (JySx) (here s are 
Pauli matrices for spin). In addition, they also induce a 
spin-independent electrostatic potential V{t) which ac- 
counts for extra scalar scattering. However, adatoms on 
graphene often tend to cluster due to ripples (as it is the 
case of H [33]) or due to a low absorption energy [34j . 
In such cases, the intervalley scattering can be neglected 
and, in the long wavelength limit, assuming potentials 
have radial symmetry, the scatterer is described by 



Vad(r)-VsoW + T/o(r), 



(1) 



where r = |r|, r being the charge carrier position. Thus, 
for r ^ R, where R is the range of the potential Vad, the 
electronic wavefunction around the K point reads: 



V'A,k,s(r) 



Akr cos 



\s) + 



(2) 



where A = ±1 indicates the carrier polarity (electron or 
hole) with energy e = Xhvpk, the ket \s = ±) describes 
the orientation of the spin along the z-axis, perpendic- 
ular to the graphene plane (s = — s); /'"*(6') and ^"^{9) 
are the elastic and inelastic ('spin-flip') scattering ampli- 
tudes at scattered angle 6, respectively. The latter is re- 
lated to the T-matrix satisfying the Lippmann-Schwinger 
equation, T(e) — Vad + VadGo(e)T(e), where Go(e) is the 
Green's function whose Fourier transform reads Go(e) = 
{e-no + XiO+)-\ Thus, /-'(0) - (Mr(e)lps'), |k| = 
k = IpI, and 9 is such that cos(6') = p • k. Let us denote 
as J-^(k, p) the 4x4 matrix whose elements are /'"' (9) in 
the spin and valley subspace. We can use symmetries of 
the Hamiltonian H = Ho + Vad to constraint the general 
form of the matrix J-^(k, p). In general, the 4x4 matrix 
J^(k,p) is a linear combination of the 16 matrices SaTp 
where a, /3 = 0, x, y, z (where a = corresponds to the 
unit matrix and a — x,y,z to a Pauli matrix). How- 
ever, the assumption of no intervalley scattering implies 
that J^(k,p) commutes with t^sq. This reduces the set 



of matrices to SqT^ where f3 — 0, z. If we further impose 
invariance under rotations about the axis perpendicular 
to the graphene plane, and time-reversal, we obtain 

J"(k, p) = asoTo + b{k A p)soTz 

+ C(k A p)SzTo + dSzTz 

+ eini • stq + e2(k A p) n2 • stq 

-f /i(kAp)ni -ST^ +/2n2 -8X2. (3) 

where the 2D vectors rii = k -I- p and n2 = k — p 
form an orthogonal set and k A p = sm9. The func- 
tions a,b,c,d,ei, fi [i = 1,2) are complex-valued func- 
tions of k and ^ = k • p = cos 6*. Note that e.g., for 
scatterers with intrinsic SOI, we can further require that 
the component of the spin perpendicular to the graphene 
plane (i.e., s • z = Sz) is conserved, which implies that 
Gj = fi = ^- Thus, the terms proportional to e^ and 
fi describe the spin-flip scattering. The symmetry con- 
straints reflected in Eq. (l3| can be used to show that, 
for an unpolarized incident beam of electrons, which is 



described by the density matrix pii, 



fSoTo|k)(k|, the 



spin-flip components of the scattering amplitude do not 
contribute to the skew scattering cross section of a single 
scatterer. This result is also obtained for an ensemble 
of the model scatterer studied below, for which electron 
transport is described by the Boltzmann equation. The 
collision integral of the latter is determined by the matrix 
of elements of J^(k, p). 

We now briefly explain how the spin Hall effect is 
enhanced by a single scatterer through skew scattering 
mechanism, and the important role played by resonant 
impurities in graphene, as well as the main differences 
with a 2DEG. To this end, let us consider a scattering 
center inducing (locally) an intrinsic SOI i.e., A/(r) ^ 0. 
As remarked above, this type of SOI conserves s^ and 
therefore /'"*(6') — 0. The details of the calculation of 
/'**(0) and the spin Hall angle are provided in the Sup- 
plementary Material (SM). Here it is sufficient to realize 
that owning to the structure of the spin-orbit coupling 
term A/r^cr^s^ [(VVo(r) x p) -s in a 2DEG], SOI scatter- 
ers break micro-reversibility |/''*(0)| 7^ |/**(— 0)|, giving 
rise to asymmetry in scattering events [c.f., terms pro- 
portional to k A p in Eq. (JSJ)]. SOI still preserves time- 
reversal symmetry, which then favors up/down spins 
to scatter symmetrically around the incident direction 
|/*'*(0)| = |J««(— ^)|^ thus explaining the formation of a 
net spin Hall current as depicted schematically in Fig. [l] 
Indeed, at the level of a single scattering event, the skew 
cross section 



r = 



27r 



d6lsine'|/'''*(6')|^ 



(4) 



is non-zero and has opposite signs for up/down spin ori- 
entations. A finite Is is the hallmark of skew scatter- 
ing. Clearly, the latter effect is absent in the first Born 
approximation, according to which the scattering ampli- 
tudes at angles ±9 coincide and hence Eq. Q is iden- 
tically zero. Moreover, we found that, contrary to the 
case of a 2DEG, a non-perturbative treatment of the SOI 
potential Vso is in general required and that, in certain 
cases, the distorted wave Born approximation, which can 
be successfully used to treat SOI in the 2DEG [26ll32]. 
fails to describe Is correctly. A few examples illustrat- 
ing the perturbative treatments and a discussion of their 
limitations in graphene are provided in the SM. 

As a measure of asymmetry in scattering events we 
adopt the so-called transport skewness; for an intrin- 
sic SOI scatterer, the latter is defined as 7 = P jo^ , 
where tr* — J d6{l — cosO) \f^'^{0)\'^ denotes the trans- 
port cross section for carrier with spin s [for Rashba SOI 
see discussion below Eq. Q]. The exact evaluation of 
7 (see SM) shows that: (i) 7 is generally non-zero for 
local SOIs of the intrinsic type, (ii) local Rashba SOIs 
induce I7I > provided that electron-hole symmetry is 
broken by a electrostatic term i.e., Vq ^ 0, and (iii) I7I 
is maximum near at resonances in a'' . To illustrate these 
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FIG. 2: Skew scattering induced by SOI impurities, (a) The 
skewness parameter 7 = I" /a" as a function of ao for a SOI 
scatterer with range R of intrinsic (Rashba) type [solid black 
line (dashed blue line)]. Here ao = VoR/{hvF) character- 
izes the strength of the local electrostatic potential caused by 
the impurity. We have taken kR = 0.15 and SOI strength 
AR/{hvF) = 0.06 (e.g., for R = 4nm these values correspond 
to hvpk — 0.01 eV and A = 20meV). Even larger values of 
7 can be found near sharper resonances (at larger Vq)- (b) 
Transport cross section versus ao for the same parameters. 
The peaks corresponds to the first resonance for Vo > 0. (c) 
Local band structure induced by SOI scatterers considered 
here. Dashed orange lines are guidelines to the eye represent- 
ing the bulk band structure of monolayer graphene. 



findings, we take disk scatterers producing an uniform 
potential inside a disk of radius R (see Fig. IT]), according 
to Vad(r) = {Vo + Vso)d{r - R), with Vso being intrin- 
sic or Rashba SOI with A//fl.(r) = A. In this work we 
take A ~ 10 meV which is consistent with ab initio cal- 
culations for metal atoms absorbed in graphene [24l [35] . 
The skewness of SOI disk scatterers in the vicinity of a 
given resonance is shown in Fig. [2J The function 7(Vo) 
is seen to follow a approximately symmetric (asymmet- 
ric) line-shape for intrinsic (Rashba) SOI with extrema. 
Note that 7 tends to zero (non-zero) as Vq — >■ for a 
Rashba (intrinsic) SOI. We also found that 7 is larger 
nearby sharp resonances, typically occurring at large Vq 
(not shown). It is well-known that already small doses 
of certain adatoms with large effective Vq values pro- 
duce resonances near at the Fermi level of graphene |30j 
that might dominate transport (see Ref. [36j for trans- 
port measurements in graphene covered with hydrogen). 
For dilute SOI disorder, the parameter 7 can therefore be 
seen as a figure of merit for the capability of generating 
net transverse spin currents via skew scattering from res- 
onant impurities. In fact, as shown in what follows, in the 
absence of other sources of impurities and at zero tem- 
perature, the spin Hall angle equals the transport skew- 
ness. Crucially, the results in Fig. [2] show that a large 
Vo is not a necessary condition to obtain large skewness. 



Although impurities such as H indeed induce giant effec- 
tive potentials Vq ~ 10^ eV (see Ref. |30j and references 
therein) , we expect impurities leading to strong SOI cou- 
plings Vso ^ 10^ meV — needed to obtain a sizable 7 — 
most likely produce effective Vq values below those found 
for resonant scatterers such as H. These impurities could 
be heavy atoms or metal nano-particles that hybridize 
with pz orbitals of graphene inducing relatively large lo- 
cal SOIs [HIMj- 

After analyzing the SHE due to a single scatterer, we 
next turn to the experimentally relevant situation of di- 
lute random ensemble of scatterers. To describe this sit- 
uation, we need to solve the Boltzmann transport equa- 
tions (BTEs), which for the steady state in the presence 
of a uniform electric field E = Ee^ read as [37| 

Vkn,(k) • (-eAkE) = ^M^,,.(k,p) ^.(p) - n,(k)] , 

p,s' 

(5) 
where Ak = sign(ek) is the carrier polarity, ^^(k) is 
the distribution function for spin s and Wss' (k, p) oc 
\r''i0)\^ (5(ek - Cp) with e = Z (k,p) is the quantum- 
mechanical scattering probability for processes with k — > 
p and s — > s' . Notice that, skew scattering implies that 
Wss'iKp) ^ W^ss'(p,k) [c.f., Eq. (§]. W^.s'(k,p) = 
^^^-^ W l^} (k, p) takes into account all sources of dis- 
order, where i? > 1 is the number of such sources. The 
above BTEs can be solved exactly, within linear response 
theory, using the ansatz for (5n<,(k) = ns(k) — n'^(k). 



0.02 F 



Sn^i^) = Vkn°(k) • [^s(k)eE -I- B,(k) (z x eE) 



(6) 



where ^"(k) is the Fermi-Dirac distribution. We deter- 
mine ^^(k) and i3s(k) by solving Eq. ^\ exactly in linear 
response. The charge and spin Hall currents follow from 
Eq. ([6| according to 3^ = ix,-\ +jx,i and j^vL = 3y,^-jy.x^ 
where js — — e^^, ^j. (57is(k)AkWk- Here, g^ — 2 stands 
for the valley degeneracy and v^ = (l/^)Vkek- At zero 
temperature, the spin Hall angle Osh = Jsh/Jx can be 
shown to assume the simple form 



7sH I T=0 



T± 



(7) 



where 



= Es',p(l-^''*'cos0)W^.v(k,p) and rji = 
J2pS,in9WssO^,p) are the transport scattering rate and 
skew scattering rate, respectively, evaluated at the Fermi 
energy ep. The spin Hall angle equals the averaged trans- 
port skewness 7 = P/a^ since a^ = '^a^F''^af^a — ''"ir^ 
and I^ = ^^ vpTiala = '''jT ' where n^ is the impurity 
areal density from source a. 

The relaxation rate has a distinct form in the spin- 
flip channel viz., t^^^^ = X]p (1 + cos^) Ws5(k,p), and 
hence, for consistency, in the calculation of the skew- 
ness of a Rashba scatterer (Fig. l2|, we have defined 
^' = Y.S' S de{l - ss' cos e)\J'"' {0)\'^. We believe 
that the commonly employed definition, r,, ^|^''pp''°'^-.' — 
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FIG. 3: Spin Hall angle as a function of Fermi energy for a 
dilute random distribution of intrinsic SOI scatterers. Panel 
a) shows 9sH at zero temperature for impurities producing a 
(non-disordered) local electrostatic potential Vq. Panels b) 
and c) show S^h at difFerente temperatures and considering a 
random Vb potential with uniform distribution Vo G [0, Al/]. 
In all panels we have taken A/ — 25 meV. 



J2r, (1 ^ COS 6) Wss(k,p), is a spurious result of approx- 
imate treatments of transport equations; the relaxation 
rates defined here, on the other hand, result from the 
exact solution of the BTEs (in first order in E). We 
also note that, for the type of potentials considered here, 
Wss oc |/ss(^)P is a function of cos{9), and thus the spin- 
flip amplitude does not contribute to r,^^. Further details 
of the solution and careful comparison with approximate 
treatments of the BTE will be reported elsewhere. 

Large SHEs are expected in relatively clean samples 
and when cross sections for unintentional sources of dis- 
order are off resonance. In Fig. 3 [panels a), b) and c)] we 
show the spin Hall angle [Eq. (7|] as a function of Fermi 
energy for pristine graphene decorated with a small den- 
sity of intrinsic SOIs (6'sh induced by Rashba scatterers 
are of the same order of magnitude and hence are not 
shown). The values obtained are comparable with those 
found in metals \0sh\ ^ 0.01 - 0.1 [3 H [12] and are ro- 
bust with respect to thermal fluctuations and disorder 
averaging [compare curves in a) and c)]. Spin Hall an- 
gles of the order of 0.1 at room temperature are obtained 
for 'fat scatterers' with R of the order of a few nanome- 
ters [see panel b)]. We thus expect clusters of metals or 
nano-particles deposited in graphene to drive giant spin 
Hall currents. As a consistency check, we verified that 
the presence of a localized magnetic moment at the im- 
purity does not supress the SHE when realistic values of 
the local exchange coupling are considered (see SM). 

Conclusions. We have shown that fat scatterers — 
adatoms, molecules or nano-particles — inducing effective 
local spin-orbit couplings of intrinsic or Rashba type in 
a graphene sheet originate robust SHEs via assymetric 
(skew) scattering, with spin Hall angles comparable to 



state-of-the-art values for metals. The SHEs studied here 
are robust with respect to thermal fluctuations and disor- 
der averaging. These findings suggest that functionalized 
graphene systems can be used to design spintronic in- 
tegrated circuits with SHE-based spin-polarized current 
activation and control. 
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Supplementary Material 



I. INTRINSIC SCATTERERS 

The components of the graphene spinor ^*(r) — 
{tp^{r),^g{r))'^ are decomposed in terms of radial har- 
monics: VaIi") = E™=_oo5m,±(^)e*'"^ and ^^{r) = 



T 



,5^±(r)e^("+i)« where 61k =arg(fc^ 



iky) and 



kx ± iky = e^^^{—idr ± ^dg), where m = 0, ±1,±2, ... is 
the angular momentum quantum number, ± represents 
the spin projection and A{B) denote the sublattice. 
After separating the variables of the Hamiltonian T-Iq + 

Vad = VF{Tz<yxPx + CTyPy) + G{T - R)[Vo + A/T^Cr^S^], we 

obtain two coupled first order equations for the radial 
functions g^_|_(r) and gm±{f)- Here, (t,t and s, are 
Pauli matrices acting on the sublattice, valley and spin 
degrees of freedom, respectively. Note that the spin-orbit 
interaction does not connect the states at different Dirac 
cones Tx = ±1 and therefore each valley can be treated 
separately. Using the partial-waves expansion for spin- 
dependent scattering, the asymptotic form for the spinor 
wave function for an incoming wave-function with spin s 
at Dirac point K is given by 



V'A,k,s(r) = 



gifcrcos(e)|g. 



+ 



\J—ir 



1 



e'^^ls) 



(1) 



where A = ±1 denotes the carrier polarity, s = —s and 
j^^{6) and /*''(6') are the the scattering amplitudes with- 
out and with spin- flip, respectively. 

Inside the disk {r < R), the eigenvalues of the Hamil- 
tonian arc given hy e = E — Vq — Xy^H^ypk"^ + Aj. Note 
that the spectrum has a gap of 2|A/|. The problem is 



decoupled in spin and the two-component electron wave- 
function at valley K is 



i'. 



1 



A,k.^ 



V2max(|e|,|A,|) V vVe - sAje'^- J 

(2) 
where rj — sgn(e + |A/|). We can write the spinors for the 
region inside and outside the potential as a superposition 
of angular harmonics. For electrons in the region outside 
the disk (r > R), E = Xhvpk and the partial wave m 
reads as 



'ip\.,Tn,sir,9) 



b„ 



Jm{kr)e'"'0 
/ H^\kr)e''^^ 



(3) 



whereas for r < R we have 



i>\,m,sir,0) = C„ 



imO 



y/e + sA jJ^{/3r)e 
ir]y/e - sAJ„+i(/3r)e'(™+i)« 



(4) 

where (3 — y/e^ — A'^ /{fiv-p). The boundary condition at 
r = R leads to two equations for each spin. The solution 
of these equations allow us to compute the amplitudes 
5^. After some algebra, we obtain 



Ve + sA/J,„+i(fci?)J^(/3i?) 



l^- 



sAiJ,n+l{PR)Jm{kR) 



r(i)/ 



Ve + sAiHl^l^{kR)J^{m - IV^^JAjlm+iimH^'ikR) 

I 



(5) 



The states with |e| < A/ correspond to damping waves 
since /3 is pure imaginary. It is straighfoward to show 
that the TTi-th phase-shift Sm |3] relates to the amplitude 
S^ as defined in Eq. (JSJ) according to S^ = ie*'''" sin((5f„). 

The Dirac equation for graphene in the presence of 
a scalar potential has a symmetry with respect to the 
exchange of g:^ by g?-m~i which is responsible for the 
absence of backscattering. For "H/, the symmetry is 
given by the exchange of g^ j. by g^^_i t^, which cor- 
responds to a relation (5_„j = (^m+i between the spin- 
dependent phase-shifts. The to = and to = — 1 chan- 
nels provide the main contributions for the spin depen- 
dent differential cross sections. If we consider the rela- 
tion (5_^ = (5^_|_]^ , we can write the scattering ampli- 
tudes in terms of the to = channels for both spins: 



J" [9) = v'2/7rfc(e*^o sinJ^ -t- e^^« 



sinJn 



The spin 



polarization is given by P{e) = \I+ {9)\^ - \r [9)\^ = 
— (S/ttA:) sin(i5|j" — 5q) sin(0), which is maximum for 6 — 
±7r/2. These calculations can be performed in the other 
valley where t^ = — 1 . In fact, it is straightforward 



to show that the charge/spin conductivity tensor is the 
same in both valleys. Notice that taking into account 
both valleys of graphene, the Hamiltonian reads 



+ {Vo + sAiT,a,)6{r-R), 



(6) 



where r^ — ±1 is the valley index and the above 
Hamiltonian acts on the space of spinorial functions in 

the form {r\iP) = f^f (r), *f (r), ^^'(r), *f' (r)") (^ 

\Xs)- The effective Hamiltonian admits an even sim- 
pler representation if one swaps the K' components of 
the spinorial wavefunction. Indeed, defining {r\ip) — 

f^f (r),^f (r),*f' (r),*5"(r)) (g) \xs), we obtain an 
effective representation in the form of two copies of the 
Dirac Hamiltonian Hq = HvfTz<t ■ p, plus a potential 
term with the same form in both valleys 



V=iVo + sAia,)0{r-R). 



(7) 



The differential cross section is therefore the same for 
states in K and K' Dirac points. 



II. TIME-REVERSAL SYMMETRY BREAKING 
SCATTERERS 

Time-reversal symmetry is broken by local exchange 
fields. If we take an extra term in the Hamiltonian of 
the form s/S.Bd{r — R) where As is the local exchange 
field strength. Using the basis discussed in the previous 
section, with A and B sublattices swapped in the K' 
valley, we obtain 



n = T^VF {(TxPx + CfyPy) 

whose eigenvalues (inside the disk) are 
es = E-Vo- sAb 



Aa/A, 



lT?v^k?. 



(8) 



(9) 



For electrons, the partial wave expansion for this problem 
is given by 



r^Ar,e) = C„ 



Ve, + gAj J^(/3,r)e""^ 

(10) 
for R < r, and where Ct = sign(es — |A/|). The mo- 
mentum inside the disk depends on the spin orientation 

according to hvpPs = \ (-1 t — Aj. For r > R, the wave- 
function reads as 



rmA^,0) 



QST 



J™(fcr)e™« 
ir,J„+i(fcr)e*(™+i)'' 

Hi^'(fcr)e™^ 

1) 

1+1 



iT,i/i'|i(fcr)e'(™+i)» 



We verified that in general the exchange field does not 
destroy SHE even if As becomes comparable to other 
energy scales. 



III. RASHBA SCATTERERS 

If we consider a scatterer with a Rashba interaction, 
the diagonalization of the Hamiltonian H inside the 
disk {r < R) yields to eigenvalues e^_s = E^^g — V — 
^Afl -|- Xyjh'^vl^k'^ + A|, where ^ = ± corresponds to 
the pseudospin degeneracy [T] . For simplicity we restrict 
the subsequent analysis to carriers with positive polarity 
(electrons), A = 1, and assume |e| > 2/S.ii. The eigen- 
states of % for the K cone can be written as 



V-^, s(r) = 



1 



hvfpk 



9. It>+*e 



hvpk 



i> 



(12) 



Differently from the intrinsic SOI, Rashba SOI entangles 
spin and and pseudo spin, which means that Sz is not 
a good quantum number and spin-flips must be taken 
into account when calculating the scattering amplitudes 
(/*'* (0) ^ 0). In the case of scattering involving spin- 
fiips, it is more convenient to work directly with partial- 
wave amplitudes. As before, we can write the eigenstates 
of the electrons for the region inside and outside the 
potential as a superposition of angular harmonics. For 
r > R and a spin-up incident beam, the incoming and 
outgoing waves can be written as 



r^Ar^e) 



iJ„,+i(fcr)e'(™+i) 



t) 



^" UH(:i"(fcr)e'(™+i)« j'^^ 



S. 






For r < R, the wave function regular at the origin can 
be written as the superposition 



rmJr.e) 



i 



Q 



^rn 






''■hVF0i^'^™+'^ 



,i{m+l)e 



t> 






;> 



(14) 



where (3^ — \/e{e — 2S,lS.ji)/hvp. The matching condi- 
tions at r = i? give us four equations 



J,„(fci?) -f S}^H^^\kR) = J2 CimJmil^^R), 

(11) JrnikR) + SJ^H^^^kR) = Yl Ci^JUkR)^ 



(15) 



(16) 



eCf 



^;^ff(i)(fci?) = ^e^^™+i(/35i?), (17) 

SV[H';^X2{kR) = Y, eQm Jrn+2(/3ci?). (18) 



We remark that Eq. ( 13 ) obeys the required boundary 



conditions. By taking a superposition of partial waves 
^ ~ Sm *'"V'm we obtain the asymptotic limit of the 
scattering solution of the Dirac equation in 2D: 



V-L.lr) 



1 I g^fcrcos^i^^ 



E ^" 



m— — oo 



e^^'Hs) 



OO 



7 . ^n 



ss ima I ^ \ p^kr\ . 

7n— — OO 



(19) 



where we can identify the scattering amplitudes in the 
elastic and spin-flip channel as 



(20) 



respectively. The above expressions allow us to calculate 
the cross-sections used in the Boltzmann transport equa- 
tion, namely the transport cross section and the skew 
cross section, 




at:,^,^ = J d9{l-ss' cos 9)\r^ {9)\' , (21) 
a:i^ ^ fd0 sine\r'{e)\\ (22) 



respectively. Inserting Eqs. (pl))-(|22|) into the above ex- 
in the u! 

oo 



pressions, we obtain the useful relations: 



?iss |2 



transp 



skew 



k 



ss'Re[S:,t (5-+i)*] ,(23) 



?n— — OO 



- J2 Im[5-'(5-;,)*]. 



(24) 



m— — OO 



IV. DISTORTED WAVE BORN APPROXIMATION 
AND SHE IN GRAPHENE 
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FIG. 1: Comparison between the exact value of ASm = 
(^Vo-i-so — Svo and DWBA approximation. The parameters 
are A/s^ = 0.007 eV, Vo = 0.07 eV and R = 1 nm. 

In this section, we provide a few examples of the limita- 
tions of the distorted wave Born approximation (DWBA) 
when applied to study spin Hall effect in graphene. We 
first derive the DWBA for a general class of potentials of 
the form 



Va = Va{r) + Wair)a,, 



(25) 



where Wa{r) denotes the sublattice symmetry breaking 
term. In analogy to the derivation for a scalar potential 
in the Schrodinger equation [2], it is necessary to write 
two copies of the Dirac equation corresponding to differ- 
ent potentials, Vi and V2. The scattering amplitudes (or 
phase-shifts, 6m) of the simpler problem ?^i = ?^o + Vi 
are assumed to be known. Let us denote the eigenstates 
of Ha = Ho + Va by 



*™(^,0) 



imcfi 



Fair) 
Gair)e'' 



(26) 



and a — 1,2. We aim at finding the phase-shifts in- 
duced by the sublattice breaking term W{r). Inserting 
the ansatz ( 26 1 into the Dirac equation, and using the 



asymptotic form of graphene wavefunctions 



^mir,^) 



2 

Trkr 



cos (kr - Xm + 5"„) 
Aisin(fcr - Xm + Sm) 



(27) 



we find 
2A 



^sin(^W-€l 



drr 



/o 

hvp 



SV+ 



Fi{r)F2{r) 



Gi(r)G2(r) 



(28) 



where 5V± = V2 ±W2 ~ [Vi ±Wi) . The above equation 
is still exact. The DWBA is derived by using the "Born 
approximation" F2{r) ~ Fi{r) and G2(r) ~ Gi{r). Spe- 
cializing to the case of interest, i.e., Vi = V2 = Voif), 
Wi = 0, and W2 = W{r), the DWBA reads as 



AJ„ 



Xnk 
2hvF 



drrWir)[fmir)^+gmir)^], (29) 



?[i] 



<:m 



?[!] 



is the correction to the m- 



where A5„ 

th phase-shift Sm introduced by the sublattice breaking 
term W{r). In the above, fm{gm) are the partial-wave 
amplitudes of the simpler problem Hi = Ho + Vo{r), 
and A — sgn{E — Vg). We can use the equation above 
to calculate A6m explicitly for a disk SOI scatterer with 
Voir) = Vo0{r - R) and Wir) = Ais^9{r - R). In this 
case 



AS„ 



Xnksz Ajk 



aR 



duu [Jm{u) 



Jm+liu)^] 



2 hvpOL^ 

(30) 
where a — \k — V/hvpl- The above expression can be 
further simplified using the properties of Bessel functions 
(not shown). In Fig. [I] we can see the comparison between 
this approximation and the exact result. 

The skew cross section can be easily calculated under 
the DWBA approximation: 



„DWBA ^ Y^ /AX a;c 



771 + 



i)cos[2(ei-C+i)] 



(31) 



9 



The DWBA seems promising to compute phase-shifts 
for intrinsic SOI scatterer in the presence of a scalar 
potential However, it fails to correctly describe the 
skew cross section (and thus SHE) for other symme- 
tries or in the presence of resonant scattering. Here, we 
briefly discuss a few situations where the approximation 
is not valid. Our first example is provided by a void in 
graphene, which is described by the boundary condition 
requiring that the A-sublattice component of the spinor 
'0>, k s(r) vanishes at r = _R. Hence, the spin-independent 

part of the scattering phase shift Sin fullfils: 



tan (5,111 = 



?[!] 



JnAkR) 
Y„,{kR)- 



(32) 



Note the symmetry S\tln — ^m ■ If we assume that 
the SOI potential only acts in the edge of the void 
i.e., W{r) = RAi6{r-R)T;,a^s^, then, the DWBA yields: 



A5„ 



rdr¥^{r,,p)W{r)i;,ni'^,(t>) 



oc A/, 



[Y„,{kR)J„,+i{kR) - Y„,+i{kR)J„AkR)Y 



Y^{kR) 



(33) 



Bessel function, which implies that Ym{x)Jm+i{x) — 
Y„i+i{x)Jm{.x) = 2/{-Kx). Hence, within the DWBA, 

[21 [2] 

5m = 5-m, that is, the same symmetry as for the void 
potential, which implies the absence of skew scattering 
and therefore SHE. 



A second example is provided by a generic Rashba scat- 
terer, for which W{r) = Afl;(r) [r'^a^sy — a^s^). It can 
be shown that within the DWBA, and at the lowest or- 
der in Afl, only the spin-flip amplitude /**(0) ^ gets 
corrected and therefore the skew cross section (and hence 
SHE) is zero in this approximation, just as in the previ- 
ous example. 
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